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SEMIPARAMETRIC ESTIMATION OF FRACTIONAL 
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By Willa W. Chen 1 and Clifford M. Hurvich 

Texas A&M University and New York University 

We consider a common-components model for multivariate frac- 
tional cointegration, in which the s > 1 components have different 
memory parameters. The cointegrating rank may exceed 1. We de- 
compose the true cointegrating vectors into orthogonal fractional 
cointegrating subspaces such that vectors from distinct subspaces 
yield cointegrating errors with distinct memory parameters. We esti- 
mate each cointegrating subspace separately, using appropriate sets 
of eigenvectors of an averaged periodogram matrix of tapered, differ- 
enced observations, based on the first m Fourier frequencies, with m 
fixed. The angle between the true and estimated cointegrating sub- 
spaces is o p (l). We use the cointegrating residuals corresponding to 
an estimated cointegrating vector to obtain a consistent and asymp- 
totically normal estimate of the memory parameter for the given 
cointegrating subspace, using a univariate Gaussian semiparametric 
estimator with a bandwidth that tends to oo more slowly than n. 
We use these estimates to test for fractional cointegration and to 
consistently identify the cointegrating subspaces. 

1. Introduction. Fractional cointegration has been the subject of much 
recent attention; see, for example, the work of Robinson [16], Robinson and 
Marinucci [19], Robinson and Marinucci [18], Chen and Hurvich [3]. All 
of these papers assume either that the observed series is bivariate or that 
the cointegrating rank is 1. Arguably the most interesting case, from an 
econometric point of view, is the situation where the series is multivariate 
and has cointegrating rank which may exceed 1. This situation was cov- 
ered by Robinson and Yajima [20], who considered methods for determining 
the cointegrating rank, and also by Chen and Hurvich [4], who focused on 
estimation of the space of cointegrating vectors. 
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Chen and Hurvich [4] studied the properties of eigenvectors of an av- 
eraged periodogram matrix of differenced, tapered observations, averaging 
over the first m Fourier frequencies, where m is held fixed as the sample 
size grows. They showed that the eigenvectors corresponding to the r small- 
est eigenvalues (where r is the cointegrating rank) lie close to the space of 
true cointegrating vectors with high probability. They also presented an em- 
pirical analysis of fractional cointegration in US interest rates for bonds of 
seven different maturities. They found evidence that the cointegrating rank 
was greater than one and, furthermore, that the memory parameter of the 
cointegrating errors may take on a variety of values that differ substantially 
if cointegrating vectors corresponding to substantially different eigenvalues 
are used. This last finding, while of apparent interest from an econometric 
point of view, could not be explained directly from the theoretical results 
presented in [4] since they did not attempt in their theory to separate the 
space of cointegrating vectors into subspaces yielding different memory pa- 
rameters. 

The goals of the present paper are to exhibit a model that allows us to 
highlight these subspaces, to show that the subspaces and their correspond- 
ing memory parameters can be estimated individually and to show how to 
use the residual-based Gaussian semiparametric estimates of the memory pa- 
rameters to consistently identify the cointegrating subspaces and to test for 
fractional cointegration. By contrast, Chen and Hurvich [4] did not consider 
either testing for cointegration or estimation of the degree of cointegration. 

We first present, in Section 2, a semiparametric common-components 
model in which the components have different memory parameters, while 
the entries of the observed multivariate series have just one common mem- 
ory parameter. Next, we show that the space of cointegrating vectors can 
be decomposed into a direct sum of orthogonal cointegrating subspaces such 
that vectors from distinct subspaces yield cointegrating errors with distinct 
memory parameters. 

In Section 5, we show that each of these cointegrating subspaces can be 
separately estimated using sets of eigenvectors of the averaged periodogram 
matrix. Since m is held fixed, we are able to obtain a rate of convergence 
for the estimated cointegrating vectors that depends only on the difference 
between the memory parameters in the given and adjacent subspaces and 
is not hampered by the rate of increase of m, as in other related work (cf. 
[19], in the bivariate case). 

To each true cointegrating subspace, there corresponds an estimated coin- 
tegrating subspace spanned by an orthonormal set of eigenvectors of the 
averaged periodogram matrix, where membership in the set is determined 
by a partitioning of the sorted observed eigenvalues into contiguous groups 
of sizes that match the dimensions of the corresponding true cointegrating 
subspaces. In Section 4, we show that the eigenvalues for the kth estimated 
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cointegrating subspace are O p {n 2dk ), where n is the sample size and dk is the 
memory parameter of the cointegrating error for the kth true cointegrating 
subspace. This result and its refinements play a key role in our subsequent 
theory. 

We will show in Theorem 1 that any vector in the feth estimated cointe- 
grating subspace is, with high probability, close to the kth true cointegrating 
subspace, in the sense that the norm of the sine of the angle between these 
two subspaces converges in probability to zero. The norm of the sine of this 
angle is O p (n~ ak ), where is the shortest distance between the memory 
parameters corresponding to the given and adjacent subspaces. This implies 
that the sine of the angle between any vector in the kth estimated cointe- 
grating subspace and the kth true cointegrating subspace is O p (n~ ak ). (We 
provide more details on the notion of the sine of the angle between subspaces, 
and also the sine of the angle between a vector and a subspace, in Section 5.) 
This convergence rate, which improves as at increases, is at least as fast as 
the rates obtained for existing semiparametric estimators of cointegrating 
vectors in the bivariate case (see, e.g., [19] and the discussion in [3]), but 
not as fast as the parametric rate obtained by Hualde and Robinson [7] of 
O p (n -1 / 2 ) in the bivariate asymptotically stationary case if the difference 
{oil) between the memory parameters of the observed series and the coin- 
tegrating error is less than 1/2. Furthermore, we show in Lemma 15 that 
the normalized eigenvectors of the averaged periodogram matrix converge 
in distribution to random vectors that lie in the corresponding cointegrating 
subspace. 

We then show in Section 6 that the cointegrating residuals corresponding 
to an estimated cointegrating vector can be used to obtain a consistent and 
asymptotically normal estimate of the memory parameter for the given coin- 
tegrating subspace, using a univariate Gaussian semiparametric estimator 
with a bandwidth that tends to oo more slowly than n. We also describe a 
procedure for consistently identifying the cointegrating subspaces, that is, 
for determining the number of subspaces and their dimensions. In Section 7, 
we provide a test for fractional cointegration which is appropriate for our 
model. 

2. A fractional common components model. Suppose that the original 
data are a q x 1 time series such that the (p — l)th differences {yt} are 
weakly stationary with a common memory parameter cLq G (— p-\- 1/2, 1/2), 
where p > 1 is a fixed integer. The use of (p— l)th differences converts any 
additive polynomial trend of order p — 1 in the original series into an additive 
constant. The value of this constant is irrelevant for our purposes since the 
estimators considered here are functions of the discrete Fourier transform 
at nonzero Fourier frequencies. We can, therefore, take the mean of {yt} 
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to be zero, without loss of generality, and our estimators are invariant to 
polynomial trends of order p — 1 in the original series. 

In order to guarantee that the cointegrating relationships in the stochastic 
component of the levels are preserved in the differences, we apply a taper to 
the differences, that is, we multiply the differences by a sequence of constants 
prior to Fourier transformation. This prevents detrimental leakage effects 
due to potential overdifferencing and allows us to obtain uniform results 
over a wide range of memory parameters. A convenient family of tapers for 
use on the differences, and which we will use here, was given in Hurvich and 
Chen [8]. The exact form of the taper is given below. 

The fractional common-components model for the (q x 1) series {yt} with 
cointegrating rank r (1 < r < q) and s cointegrating subspaces (1 < s < r) is 
given by 

(1) yt = Aor4 0) + Aiu^ + • • • + A s u { t s) , 

where A k (0 < k < s) are q x a k full-rank matrices with ao = q — r and 

a\ H \-a s = r such that all columns of Ao, . . . , A s are linearly independent, 

and }j k = 0, . . . , s, are (a^ x I) processes with memory parameters 
{dfc}| =0 with — p + 1/2 < d s <••■< do < 1/2. Equation (1) can be written 

as 

(2) yt = Az t , 

where z% = vec(u|°^ , . . . , u^) and A =[Ao • . . A s ] . We will make additional 
assumptions on {z t } in Section 3. These assumptions guarantee that {z t } is 
not cointegrated. The methodology presented in this paper does not require 
either r or s to be known. 

Remark 1. Our assumption that all entries of {yt} have memory pa- 
rameter d implies that all rows of Aq are nonzero. The model (1), without 
the assumption that all entries of {yt} have a common memory parameter, 
could also be entertained (though we do not pursue this here) and would 
then include the model considered by Robinson and Yajima [20]. 

Next, we exhibit the cointegrating subspaces. For any matrix A, let A4 (A) 
denote the column space of A and let .M^(A) denote the orthogonal com- 
plement of A. Note that for k = 1, . . . , s, 

M ± (A ,...,A fc )cM ± (A ,...,A fc _ 1 ). 

Let £>o = A4(Aq) and k = 1, . . . , s, be the subspace such that 

M ± (A Q , A fc _i) = M ± (A , ...,A k )®B k 
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and B k -L M. (Ao, . . . , A k ). Hence, a nonzero vector /3 G G {1, ... , s}, 

satisfies /3'A^ = 0, I = 0, . . . , jfe - 1, and /?A fc / 0. Also, Bj _L £ fc for j + k, 
(j, k) e {0, . . . , s} and 

(3) R f/ = 5 

It can be seen from (1) and the preceding discussion that any nonzero vector 
j3 € Bk with fe € {1, . . . , s} produces a cointegrating error series {(3'yt} with 
memory parameter d k . Thus, B\,...,B S are the cointegrating subspaces. 
The space Bq, on the other hand, is the space spanned by any basis of non- 
cointegrating vectors in M 9 . Equation (3) shows that M q may be written 
as a direct sum of the space of non-cointegrating vectors and the space of 
cointegrating vectors, and that the latter space may be further decomposed 
into a direct sum of cointegrating subspaces. 

3. Assumptions. Here, we specify a linear model for the series Zt = 
vec(u[°\ . . . ,u\ ). As stated in the previous section, we assume that {u^}, 
k = 0, . .. ,s, are {a k x 1) processes with memory parameters {d k } s k=Q with 
-p + 1/2 < d s < ■ ■ ■ < d < 1/2. Define iV = {1, . . . ,a } and N k = {(a + 

• — h afc_i) + 1, . . . , (ao H h a k )} for k = 1, . . . , s. Our results in this paper 

assume s > 0, unless explicitly stated otherwise. 

Let ip k be a sequence of q x q matrices such that 

i> k = ^j\ iku3 <S>{uj)duj, 

where for each uj G [— vr,7r], *&(uj) is a complex- valued matrix such that 
S&(— uj) = &(uj) and where ipo is an identity matrix. 
Define the q x 1 vector process {zt} as 

oo 

(4) 3* = ^ ipk£t-k, 

k=—oo 

where {et = (st,i, ■ ■ ■ ,£t,q)'} ~ iid(0, 27rS), S is a symmetric positive definite 
matrix with entries a a b, a, 6 G {1, ... , q}, and -E||et|| 4 < oo, where || • || denotes 
the Euclidean norm. The spectral density matrix of {zt} is 

f(w) = ¥(u/)E¥*(c«;), cjG[-vr,vr], 

where the superscript * denotes conjugate transposition. We further assume 
that for uj G [— vr,7r], the (a,6)th entry of St?(uj) is given by 

(5) = (1 - e- i »)- d "H ah {ujyt>^\ 

where d aa = d k for a G N k , d ab < mm(d k ,d h ) for a G iV fe , 6 G N h , b / a, 
k,h = 0, . . . ,s, and for all a, 6 G {1, ... , g}, r a j,(-) are positive even real-valued 
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functions and <j) a b{') are odd real-valued functions, all continuously differ- 
entiable in an interval containing zero. It follows from (5) that the first 
derivatives of ^ a b(uj) satisfy 

(6) ^)=0(\^aa(^ bb H\ 1 / 2 \u;\- 1 ). 

In keeping with (5), we assume that we can write the spectral density matrix 
of {zt} as 

(7) f(w) = T(w)ft(w)T* 



where T(w) = diag{(l - e^) - *, . . . , (1 - e ~ iuj )- d ° , . . . , (1 - e~^)- da , . . . , (1 - 
e~ lu ')~ ds }, that is, the ath diagonal entry is (1 — e~ tu) )~ dk for all a G N k , 
k = 0, . . . , s, and 



(8) ft(w) = *t*( w )s^t(, 



is positive definite, Hermitian, continuous at zero frequency and, therefore, 
real-valued at zero frequency with ^ a u(w) = T a i ) (uj)e' t ^ ab ^ . Thus, {zt} is not 
fractionally cointegrated (see [18]). 

4. The averaged periodogram matrix and its eigenvalues. For any vector 
sequence of observations {£t}" =1 , define the tapered discrete Fourier trans- 
form by 



where Uj = 2irj/n is the jth Fourier frequency and {/ij} is the complex- 
valued taper of Hurvich and Chen [8] , 

h t = 0.5(1 -e i27Tt/n ), t = l,...,n. 

Note that p = 1 yields the no-tapering case. Next, define the tapered cross- 
periodogram matrix of two vector sequences {£,t}t=i an d {CtltLi by 

I^(u j ) = J^u j )J^(u j ). 

We will work with the (real part of the) averaged periodogram matrix of a 
sample of n observations {yt}t=i> 

m 

I m = J2Re{lyy(uJ j )}, 

where m is a fixed positive integer, m > q + 3. (This condition is motivated 
in the proof of Lemma 8.) 
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Define I m (£,t,(t) = X)j=i ^ e {-^c( w j)}- We ^ rs ^ f° cus on the asymptotic 
distribution of I m (zt,zt). Define the function (for 16E) 

A p (,)=( 2 ^ 1 2 )~ 1/2 E(^ 1 )(-l) fe A(x + 2 7 rfc) ; 
v y J k=o v y 



where 



Now, define 



A(x) -- 



2tt ix 



Vj(x) = -[A p (-x + 2nj) + A p (x + 2vrj)], 



Uj(x) = ~[A p (-x + 27Tj) - \{x + 2vrj)]. 

Define the Hermitian positive definite qx q matrix- valued measure Go on 
R by 

(9) G (dx) = n(x)f t (0)n*(x) dx 



for x > and Go(— dx) = Go(dx), where 

IL(x) = diag(e- i7rdo/2 |x|- do , . . . , e- i7Tdo/2 \x\- do , . . . , e - i7rd ° /2 \x\- d ° , . . . , 

e -ind B /2\ x \-d s \ 

Let U n and V n he qx m matrices given by 

(10) U n = d~ 1 Re(J Zil ,...,J Z)m ) and V n = d" 1 lm(J Zjl , . . . , J z , m ). 

Lemma 1. Let d n be a (qx q) diagonal matrix with ith diagonal entry 
n dk , i£N k , k = 0,...,s and Q n = d' 1 I m (z u Zt)d~ l = (U„, VJ(U n , V n )'. 
If m> q, then 

Qn-^UU' + W, 

where U =(JJ\, . . . , U m ) and V ={V\, . . . , V m ), Uj, V k are qx 1 vectors and 
vec(U,V) is a 2mq-variate normal random variable with zero mean and 
covariance matrix S determined by 



E(U j U' k )= / V j(x)v k (x)G (dx) 



nVjVl) = / vAx)v k {x)G Q {dx), 



nUjVi) = / v 3 (x)v k (x)G (dx). 



Furthermore, UU' + W is positive definite and has distinct eigenvalues 
with probability 1. 
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Proof. The proof is identical to the proof of Lemma 1, Corollary 1 and 
2 of [4]. □ 

We next derive upper and lower bounds for the eigenvalues of I m (yt,yt)- 
We will use the notation Aj(-) for the jth eigenvalue of a given Hermitian 
matrix, Xj(-) > Aj+i(-). Also, we let Xj = Xj(I m (yt,yt))- We have the follow- 
ing lemma: 

Lemma 2. Xj = O p (n 2dk ) for j eN k , k = 0,...,s. 

In the case k > 1, the upper bound in Lemma 2 strengthens Lemma 4 
of [4]. 

Lemma 3. Let jt = max{j : j £ N k } and let Qn^ be the leading j£ x j| 
principal submatrix of Q n for k = 0, . . . , s. Then 

n-^X 3l >c k X 3l (^)^^, 

(k) 

where c k > and rn* is a random variable that has no mass at 0. 

•* k 

5. Estimation of the cointegrating subspaces. Let X(-) = [xi(0 ■ ■ ■ Xq{')] 
be an orthogonal matrix such that Xj(') is the eigenvector corresponding 
to the jth largest eigenvalue Xj(-) of a given symmetric q x q matrix and 
let Xfc(-) be a matrix with columns 3 e N k , f° r k = 0, . . . , s. Also, 

let Xj = Xj(Im{yt,yt)), X = X(I m (y t ,y t )) and X k = X k (I m (y t , y t )). For k = 
0,1, ... ,s, let B k be a g x a k matrix with orthonormal columns such that 
M{B k ) = B k and let B = [B ... B s ]. Since B'B = I, it follows that for 
any q x q matrix P, B'PB is similar to P, that is, A,(P) = A,(B'PB) and 
Xj (P)=B' Xi (B'PB). 

Define 

& = B'l m (y t ,y t )B 

and partition 3> into (s + l) 2 blocks such that the (k,£) block <f> k g has 
dimension (a k x a/) for k,t = 0, . . . , s. Define <&d = diag[<&oo> ■••■> *&ss] and 
A* = 3> - so that 

$ = <& D + A*. 

We have 

I m (yt,yt) = B*B' = B$ D B' + BA*B' =: H+AH, 

so we can think of I m (yt,yt) as a perturbed version of H. Using results of 
Barlow and Slapnicar [2] on perturbation theory for eigenvalues and eigen- 
vectors of nonrandom Hermitian matrices, we will show in Lemma 4 that 



FRACTIONAL COINTEGRATING SUBSPACES 



9 



the &;th estimated cointegrating subspace A4(X. k ) is close to M(X. k (H)) in 
the sense that the norm of the sine of the angle between the two subspaces 
converges to in probability. 

Let &(•,•) denote the matrix of canonical angles between two subspaces 
of the same dimension (see, e.g., [22], page 43). The notion of the sine of 
the angle between two subspaces of the same dimension is given in [5]. 
For simplicity, suppose that S and T are both real gxa matrices (q > a) 
with orthonormal columns. Then the orthogonal projector into M(T) is 
given by TT' and the projector into the orthogonal complement .M^(T) 
of .M(T) is given by I — TT', where I is a q x q identity matrix. The 
sine of the angle between .M(S) and M(T) is an a x a matrix defined 
in [5] and denoted by sin6(.M(S), M(T)). It follows from [5], page 10 that 
\\smQ(M(S),M(T))\\ F = ||(I - TT^SS'I^, where \\-\\ F is the Frobenius 
norm. It follows from [22], Corollary 5.4, page 43 that 



where T is a matrix with orthonormal columns spanning A4 (T), so that 
|| (T J -) / S is the square root of the sum of the squared lengths of the residu- 
als from the orthogonal projections of the columns of S on the space M(T). 

For any nonzero vector x 6 M(S), the sine of the angle between x and 
the subspace M. (T) is a real number defined as 



In Lemma 5, we show that under the additional assumption that the pro- 
cess is Gaussian, .M(Xfc(H)) is equal to B k with probability approaching 
one, for k = 0, . . . , s. Lemmas 4 and 5, taken together, imply our Theorem 1, 
stating that if the process is Gaussian, then the fcth estimated cointegrat- 
ing subspace Ai(X. k ) is close to the corresponding true cointegrating sub- 
space B k , in the sense that || sm@{Ai(K k ),B k }\\F = O p (n~ ak ), where a& is 
the shortest distance between the memory parameters corresponding to the 
given and adjacent subspaces, that is, 



(11) 



sinG(^(S),A^(T))|| F = ||(T ± ) , S 




II fl - TT') 

\\x\\ 



see [24], page 274. It then follows from (11) that 

max .\ S md(x,M(T))\<\\(T ± yS 

x£M{S) 



= < 



d - di, 

min{(c4_i - d k ), (d k - d k+ i)} 
d s -i — d s , 



k = 



k = 1, . . . , s — 1 
k = s. 



Lemma 4. The sine of the angle between M(X. k ) an d -M(X/%(H)) sat- 
isfies 

\\sme{M(X k (U)),M(X k )}\\ F = O p (n-^). 
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The following Gaussianity assumption is sufficient for obtaining a rate 
at which P(.M (X&(H)) / B k ) converges to zero. More specifically, the as- 
sumption allows us to bound the inverse second moment of eigenvalues of 
Q n . We believe that such bounds, and therefore Lemma 5, hold without the 
Gaussianity assumption, but we will not pursue this here. 

Assumption 1. The process {et} in (4) is Gaussian. 

Lemma 5. Under Assumption 1, P(M(X k (H)) / B k ) = 0(n~ 2oik ), k = 
0,...,s. 

The following theorem is a corollary of Lemmas 4 and 5: 
Theorem 1. Under Assumption 1, 



6. Estimation of the memory parameters using cointegrating residuals. 

Let b = Xa, where a £ {1, ... ,q}. Recall that Xa is the eigenvector of I m (yt,yt) 
corresponding to the ath largest eigenvalue of the matrix. Then there exists 
a uniquely defined value k G {0, . . . , s} such that a G N k . Note that k is fixed 
but unknown. We then use this vector b to construct the residual process 
{vt}, where 



(12) v t := b'y t = b'Aou^ + b'A lU ^ + ■■■ + b'A k u\ k) + ■■■ + b'A s u ( t 



The periodogram of {v t } is I vv (ujj) = b' AI zz (u>j)A'b. We consider the Gaus- 
sian semiparametric estimator or GSE (see [10, 17]) for d aa [see (5)] based 
on {v t }, 



where 6= [Ai, A 2 ], -p + 0.5 < A x < A 2 < 0.5, uj- = 2-irj/n, ] = j + (p- l)/2 



shifts inherent in our tapering scheme and thereby reduce finite-sample bias, 
as was also done in [8]. 

The two theorems below establish the consistency and the limiting distri- 
bution of the d aa , under some additional conditions on the transfer function 
tfjjfw) =T a &(w)e i *«»H; see (5). Following [9], we define a smoothness class 



smG{M(X k ),B k }\\ F = O p (n 



k = 0, . . . , s. 




and 
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for transfer functions as follows. For p> 1 and 1 < p < 2, let C*(p,p) be the 
set of continuously differentiable functions u on [— tt,tt] such that for all x,y 
with \x\ £ (0, 7r], \y\ 6 (0,7r], 

max < z < n \u(z)\ \u(x)-u(y)\ \y - A 

— i TTT — ^' 



min < z <^|u(z)| min < 2 < T |u(z)| min(|a:|, ' 



I it' fx) — \y — xK p ^ 



mino^z^Tr |it(z)| [min(|x|, |y|)]^' 

It follows from the discussion in [9] that if ^> aa (ui) is the transfer function 
of a stationary and invertible autoregressive moving average process, or of 
the short-memory component of a stationary and invertible fractional Gaus- 
sian noise with a suitable choice of the moving average representation, then 
\P]j (u;) G C*(p, p) for some p, with p = 2. 
We now state an assumption on ^ . 

Assumption 2. For all a, be {1, ... ,(/}, G £*(p,p) for some p > 1 
and some /? S (1, 2]. 

Note that this assumption is global in that it pertains to the behavior 
of at all frequencies. By contrast, our estimation of the d aa is based 
on frequencies in a shrinking neighborhood around zero. It seems plausible, 
then, that a local version of Assumption 2 would suffice for our purposes, 
although we do not pursue this here. 

The following standard assumption is needed to establish the consistency 
of d aa : 

Assumption 3A. As n^oo, 

1 m n 
— + — ^0. 

m n n 

Theorem 2. Under Assumptions 1, 2 and 3A, for a S {1, . . . , q}, d aa — ► 



The next assumption is used for establishing the asymptotic normality of 

m n (d a a ~ d aa ). 

Assumption 3B. Suppose that a e N k . 
(i) If k e {1, . . . , s}, then 4_ x - d k > 1/2. 
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(ii) If k € {0, . . . , s — 1}, then as n — ► oo ; 

TO n n 2(c( fc -d fe+ i) 

Note that part (i) is vacuous if k = and part (ii) is vacuous if k = s. 
Assumption 3B may be compared with the assumptions in Theorems 2 and 
4 of Velasco [23] which the author required for residual-based estimators of 
the memory parameters of a bivariate fractionally cointegrated system. The 
problem here is that a linear combination of series with slightly different 
memory parameters will typically have an irregular short-memory compo- 
nent in its spectral density. 

To present the asymptotic variance of d aa , we define 

_ T(Ap-3)T\p) 

p r 4 (2p-i) ' 

Theorem 3. Under Assumptions 1, 2 and 3B, for a £ {1, . . . ,q}, 

my 2 (4a-daa)^iV(0,V 4 )- 

1/2 * 

Note that in Theorem 3, the limiting distribution of m n (d aa — d aa ) has 
mean zero. This asymptotic unbiasedness is ensured by Assumption 3B, 
which places strong restrictions on the separation between the memory pa- 
rameters and also places a potentially stringent upper bound on the band- 
width m n . A much weaker and, indeed, more standard assumption involving 
only m n is the following: 



Assumption 3C. As oo, 

1 | m l + 2 nog 2 m n Q 
m n n 2 P 

If we account for the asymptotic bias, which can be determined from 
Lemma 20, and use Assumption 3C, we obtain the following result: 

Corollary 1. Suppose aGN^, where k G {0, . . . , s}. Under Assump- 
tions 1, 2 and 3C, we have 

ml/ 2 (d aa - ~ Vn) ~^ N(0, $ p /4), 

where fJ, n = O p (mn k dk ~ 1 Iwi dk + 1 ^ ^ e o p (m n k dk ~ 1 ) term is vacuous if 
k = and the O p {uJm n dk+1 ) term is vacuous if k = s. 
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Here, we present some results on the vector of GSE-estimated memory 
parameters, d = (du, . . . , d qq )' , which is an estimate of d = (du, . . . ,d qq )'. 
Let wt = X'yt be the q x 1 residual vector so that the entries of d are based 

on those of Wt . Note that by Lemma 15, X X, where X is a continuous 
function of U and V in Lemma 1. We will need the following assumption 
for our results: 

Assumption 3D. (i) For all fee {0, ...,s}, a k > 1/2. 
(ii) As re — ► oo, 

1 rnl+ 2 ^ log 2 m n ^ o 
m n re 2 ? 

where £ = minjminfc o^, p}. 

Corollary 2. Under Assumptions 1, 2 and 3D, 

m][ 2 {d -d)^N (o, ^(diag ft)- 1 oft oft o (diag , 

where 

fi = E(X'Af t (0)A'X). 



Remark 2. Simulation results not shown here reveal that the small- 
sample bias is reduced and the variance is stabilized if the GSE estimators 
omit the first m + p — 1 frequencies. This does not affect the validity of 
Corollary 2. Note that if no frequencies are omitted, then the first m + p—1 
frequencies are used twice: once for estimating the cointegrating vector and 
once for estimating the memory parameter. If the frequencies are omitted, 
the finite-sample approximation to the variance in [8] is quite accurate. 

Corollary 2 yields the following result on the asymptotic distribution of 

1/2 ~ 

fn-n {d aa — d-bb — [d aa — dbb)), under conditions that ensure asymptotic un- 
biasedness: 



Corollary 3. Under the assumptions of Corollary 2, for a, b £ {1, . . . , q}, 
a^b, 

m]l 2 {d aa - d bb - (daa ~ d bb )) iV^O, ^1 — j^-)) • 



Next, we modify Corollary 3 to include a bias term, thereby allowing for 
weaker assumptions. 
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Corollary 4. If a G Nk, b G Nh, for k,h£ {0, . . . , s}, then under the 
assumptions of Corollary 1, 

m}l\d aa - d bb - (d aa - d bb ) -Jl n )^N (o, Ml - ^-)) , 

where 

Jin = Opimi^- 1 + m^- 1 + u^^ 1 + J m h - dh+1 ). 

Given data from model (1), assumed to possess fractional cointegration, 
the number s > of cointegrating subspaces and their dimensions a± , . . . , a s , 
as well as the dimension ao of the non-cointegrating space, will be unknown 
in general. Here, we assume Gaussianity. Let 5* > be the minimum sep- 
aration between the memory parameters, 5* = min(do — di, . . . , d s -i — d s ), 
and assume first that 5* > 1/2. We can compare the GSE estimators djj and 
(ij+ij+i for j = l,...,q using a bandwidth m n satisfying Assumption 3D, 
part (ii), with £ = min{<5*,p}. Fix an e G (0, 1/2) and a C > 0. Then, for 
each j G {1, . . . ,q — 1}, we declare that djj — dj + ij + i / if and only if 

djj — dj + ij + i > Cm n l l 2+£ . This leads to a procedure for consistently iden- 
tifying s, ao, . . . ,a s , which can be justified by Corollary 3. A more compli- 
cated identification procedure, justified by Corollary 4, may be constructed 
if 5* < 1/2. This requires further tuning parameters which depend on 5*, 
owing to the presence of the nonstandard term fl n , which increases as the 
separation of the relevant memory parameters decreases. Unfortunately, such 
procedures are of limited practical value as 5* will be unknown in practice. 
Nevertheless, we note that lower bounds on 6* (assuming s = 1) arise implic- 
itly or explicitly in other works on semiparametric fractional cointegration. 
(See [20], Assumption D and [23], Theorems 2 and 4.) 

7. Testing for fractional cointegration. In model (1), used throughout 
the paper thus far, we have assumed that s > 1 so that cointegration exists. 
Here, we expand model (1) to include the case of no cointegration (s = 0, or 
equivalently, r = 0), that is, 

(14) yt = A ui°\ 

where Ao is q x q with linearly independent columns and all entries of u[°^ 
have memory parameter do. 

In practice, it is of interest to test for the presence of fractional cointegra- 
tion. Such a test was proposed by Marinucci and Robinson ([13], pages 236- 
237), following on from an idea originally suggested in a different context by 
Hausman [6] , using a comparison of two estimates of do , one based on a mul- 
tivariate Gaussian semiparametric estimator (see [11]) using {?/(}™ =1 with an 
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imposed restriction that all entries have the same memory parameter, and 
the other estimator based on a univariate Gaussian semiparametric estima- 
tor of do using (say) the first entry {yij} of {yt}- It seems possible to use this 
idea, together with differencing and tapering, to yield a test for fractional 
integration in the current context, although we do not pursue this here. We 
focus instead on residual-based methods in which estimated memory pa- 
rameters based on the various cointegrating residual series are compared. 
In a bivariate context, Velasco [23] has considered properties of semipara- 
metric memory parameter estimates based on cointegrating residuals under 
certain assumptions on the rate of convergence of the semiparametric esti- 
mator of the cointegrating parameters. However, the author did not present 
a test for cointegration since his assumptions ruled out the no-cointegration 
case. 

For our GSE estimators d based on cointegrating residuals, we have the 
following extensions of Corollary 2 and Corollary 3 to the no-cointegration 
case (14): 

Theorem 4. 
tion, then 

Corollary 5. Under Assumptions 1, 2 and 3C, if there is no cointe- 
gration, then for a, b E {1, ... , q}, 

^(4.-^-2.^(1-^)). 
Corollaries 4 and 5 justify a conservative hypothesis test for the null hy- 

1/2 " 

pothesis of no cointegration based on the test statistic T n = m n {d\\ — d qq ) 
whereby, for a nominal level a test, the null hypothesis is rejected in favor 
of the cointegration alternative hypothesis if and only if T n > ($ p /2) 1/ ' 2 z a / 2 - 
Here, a bandwidth m n satisfying Assumption 3C should be used. The test 
is conservative since (<3? p /2) is an upper bound for the asymptotic variance 
of T n . 

8. Proofs. 

8.1. Proofs for Section 4. 

Proof of Lemma 2. Note that 

* = B'A/ m (z^)A'B, 



Under Assumptions 1, 2 and 3C, if there is no cointegra- 
d)-^N (o, ^-{diag ft)" 1 o ft o ft o (diag ft)~A . 
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where B'A is an upper triangular block matrix. We first partition <& into 
(s + 1) x (s + 1) blocks such that the (k,£) block has dimension (a^ x ag). 

Let z[ k) = {u {k \...,uf ] ) and A< fc > = [A fc ... A s ], k = 0, 1, . . . , s. We have 
, . * M = B' fc A( fc )j ro (4* ) ,4 <) )AW'B / , forfc<£, M = 0, 1, . . . , s, 

Fix a value of k G {0, . . . , s}. Note that by Lemma 1, all the elements in the 
kth. block, are O p {n 2dk ). Now, 



jeNfc jeiV fe u---u7V s v=k 



n 2d ^ 



See, for example, Theorem 14 of Magnus and Neudecker ([12], page 211). 
We have Xj = O p (n 2dk ) for j £ N k . □ 

(k) 

Proof of Lemma 3. Following from Lemma 1, Qn converges in dis- 
tribution to a matrix that is positive definite with probability one. Since 
an eigenvalue of a matrix is a continuous function of the entries of the ma- 
trix, we conclude that Xj* (Qn ), the smallest eigenvalue of Qn , converges 
in distribution to a random variable that has no mass at zero. To prove 
n~ 2dk Xj* > CkXj*(Q,^), we construct another, similar, matrix for I m (yt,yt)- 



Let C s = Ai(A s ) and C k , k = 0, . . . , s — 1, be the subspaces such that 

M ± (A k+1 , ...,A S )= M ± (A k , . . . , A s ) C k 

and C k -L M. ± (A k , . . . , A s ). For k G {0, . . . , s}, let C k be a q x a k matrix with 
orthonormal columns such that M{C k ) = C k and C =[Co . . . CJ. By this 
construction, P = C'A is a lower triangular block matrix and PI m (zt, Zt)~P' = 

Ol m (y t ,y t )C := W is similar to I m {y u yt)- Let P«, W« and d {k) , k = 
0, 1, . . . , s, be the leading j% x principal submatrices of P, W and d n , re- 
spectively. Also, let z[ k) = {uf\. . .,u { t k) ), k = 0,1,..., s. Note that = 
p( fe )7 m (5 t (fc) ,5f ) )P( fc )'. By Corollary 2.2.1 of Anderson and Das Gupta [1], 

A^(WW) > X rk {I m (t\t ] )}X^ {k) '^ k) ) = c k X rk {I m (t\t ) )} 

and 

X j; {I m (4 k \~4 k) )} = A Jt! {d( fc )Q«d( fc )} > A jA! (d( fc ))A jA! (Q( fc ))A i: (dW) 



2d k \ .*(n( k ) 



Ai|(Q} 



Applying the Sturmian separation theorem ([15], page 64), we have 

J k 



Air = A i2 (W) > A,.(w( fc )) > c fc A i2 {J m (4 fc U (fc) )} > ^ 24 A,(Qi fc) ). □ 
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8.2. Proofs for Section 5. 

Proof of Lemma 4. Since X£(H) = [X (H) . . . X fe _i(H) X fc+1 (H) . . . 
Xg(H)], we have 

s 

\\sme{M(X k (H)),M(X k )}\\ F < ||(X^(H))*X fe || F < ll( x <(H))*X fc || F 

£=0,£^k 



O p (max n H* - *l) = O p (n" afc ), 



by Lemma 7. □ 

Proof of Lemma 5. For k = 1, . . . , s — 1, we have 
P(M(X fc (H)) = B fc ) 



p({MX fc (H)n ^ = o}n{^x,(H)n B* = o}Y 



Hence, 

P{M{X k (H))^B k ) 

= p({^x fc (H)n 6 f /o}u{MX fc (H)n B e ^0 

^ ^ «<fe-l ^ ^ l>k+l 

<p(W fe (H)n B^oj+pUitHjn B t ^0 

by Lemma 10. Similarly, 

P(M(X (H)) / Bo) = O p (n- 2d0+Ml ) 

and 

P(M(X s (U))^B s ) = O p (n- 2d ^ +2d °). 
We have thus completed the proof. □ 

We will need the following lemma for the proof of Lemma 7. First we 
write U and V defined in Lemma 1 as 



U = [U ... UX V = [V ... V; 

where and are a k x m matrices. 
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Lemma 6. Let K = diag(Bg A , . . . , B' s A s ) . Then 

d-^d" 1 K(UU' + W')K', 

where d n is defined as in Lemma 1. Furthermore, 

d-^d" 1 K(UU' + VV%K', 

where (UU' + VV')d = diag(U U^ + V V^, . . . , U S U' S + V S V^) andK(UXJ' + 
VV'jflK' is positive definite and has distinctive eigenvalues with probabil- 
ity 1. 

Proof. We write <& = K/ m (zj, zt)K + R, where R is a symmetric ma- 
trix with its (fc,f)th entry given by 

R^ = B' fe A fc / m (^\^ +1 V ( ' +1)/ B. + B;A( fc + 1 )j m (4 fc+1 \«f))A^B, 

+ B' fc A( fc+1 )/ m (zf +1) ,zf +1) )A ( ' +1) 'B^ 

for A; < £, £ = 0, 1, . . . , (s - 1), R fcs = B' fe A( fc+1 > J m (4 fc+1) , 4 S) ) A' S B S for it < s 
and R ss = 0. Thus, 

d-^d" 1 = d- 1 K/ m (z t , ^K'd" 1 + d^Rd" 1 . 

By Lemma 1, 

d^KI^^K'd- 1 K(UU' + W')K' 
and the (&,£)th entry of d~ 1 Rd~ 1 is 

n~ dk ~ de K k£ = O p (n*+ 1- * + n dfe + 1_dfe ) = o p (l) for fc < 1 

We have proved the first limiting distribution of the lemma. It follows that 
the feth diagonal block of has the limiting distribution 

n~ 2d ^ kk K kk (XJ k U' k + V fc V' fc )K' fcfc 

and <& kk is positive definite, having distinctive eigenvalues with probability 1 
by Okamoto [14]. □ 

Lemma 7. ||X*(H)X fe || F = O p {n~\ d "^\) for all £, k e {0, 1, s} with 
£±k. 

PROOF. Since ||XJ(H)X fc || F = ||X|(* D )B'BX fc (*)|| F = ||X|(* I? )X fc (*)|| J p, 
we prove this lemma by showing that 

||X|(* I? )X fc (*)|| F = O p (n-l dfc -^l). 

Let A = diag{A j5 j = l,...,q} and = {Xj,j G N k }. We define A(*r>) 
and A^(<&d) similarly for We will use the bound for the error in 
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two subspaces within the nonzero space from Theorem 4.1 of Barlow and 
Slapnicar [2] (which can be shown to apply in our context with probability 
one), that is, 

\\X e (* D )X k m\ F < re lga P (AW(* D ),A^)) ' 

where 



relgap(A w A (fe) ) = min 



It is sufficient to show that 

(16) ||A- 1 /2 ($d)x *($ d )A*X(*)A- 1 / 2 || f = O p {\) 
and 

(17) — n ( v -\ d k-<ii\\ 

1 j relgap(AW(* D ),A( fc )) " p{ h 

By Lemmas 2, 3 and 6, relgap(AW(# D ), A^)) = O p (ral dfc -*l) and n ^l dfc ~^l 
x relgap(A^) (<&£>), A( fc )) > q^., where qj. is a random variable that has no 
mass at 0. We thus have (17). We next prove (16). Note that by Lemmas 1 
and 6, 

d n *- 1 d n ^K'- 1 (UU' + VV')~ 1 K- 1 . 

Hence, d n X(*)A" 1 / 2 = O p (l) since d n * -1 d n = d n X(*)A- 1 / 2 A~ 1 / 2 X'(*)d n 
= Op(l). Similarly, A~ 1 / 2 (* D )X*(* D )d n = O p (l). We have 

||A- 1 / 2 (* D )X*(* D )A*X(*)A- 1 / 2 || F 

= ||A- 1 /2 ($D)x * ($D)dnd -i A$d -i dnX (*)A-i/2|| F 

< ||A- 1 / 2 (* D )X*(* D )d n || F ||d- 1 A*d- 1 || F ||d n X(*)A- 1 / 2 || F 

= o P (i), 

by Lemma 6. Hence, ||XJ(* D )X fc (*)|| F = O p (?i-^- d ^). □ 

We need the following two lemmas for the proof of Lemma 9: 

Lemma 8. Under Assumption 1, there exists a finite constant C not 
depending on n such that for all sufficiently large n, 

E[A 2 (Q" 1 )]<a 
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Proof. Note that 

Q n = (U n ,V n )(U n ,V n )', 
where U n and V n are defined in Eq. (10). Let 

T(W n ) = X 2 1 (Q- 1 ), 

where W n = vec(U n , V n ). By Assumption 1, W n ~ N(0,3 n ), where 3 n = 
cov(W n ) and H n — > H, the covariance matrix of vec(U, V) in Lemma 1. It 
was shown in [4] that S is positive definite. Thus, for all sufficiently large 
n, H n is invertible and Ai(S n ) — > Ai(S) > 0. 
For all sufficiently large n, 

E s JT(^ n )] = (2vr)- m9 |H n r 1 / 2 / T(x) e - X ' 3 ^ 2 dx. 

JM. 2m i 

Since x'H" 1 ^' > x'x/ Xi(S n ), we have 

e -x'E-V2 < e -xV2Ai(H„)_ 

Since Ai(H n ) — » Ai(H) > and since |3 n |- 1//2 — ► | S | x / 2 > 0, there exist 
constants C\ > and C2 > such that for all sufficiently large n, 

E Sn [T{W n )\ <cJ T(x)e- C * x ' x ' 2 dx = C, 

JR 2m <i 

a finite constant which does not depend on n. The above integral is the 
second moment of the largest eigenvalue of an inverse Wishart matrix and 
is hence bounded by a finite constant [21], in view of our assumption that 
m > q + 3. □ 

Lemma 9. Define Eke to be an event, Eke = {A flfc (*fcfc) > Ai < 
k < £ < s. Then under Assumption 1, 

P{E c M ) = 0{n- 2dk+2dl ). 

Proof. For i > k, I = 1, . . . , s, we have, by Chebyshev's inequality and 
the Cauchy-Schwarz inequality, 

P(E c M ) = P{\ ak (<f>kk)<\i($ee)} 

_ p [ - 2dh+2df n-^Ai(^) >t 
I n- d k\ ak {3> k k) 

< n -24+2d £E l/2 [A 2 {n -2d e ^)] E l/2 [A 2 ^ $ -l )] _ 

Since E[A?(n -Mt *#)] < E[trace 2 (?i~ 2 **«)] < C by Assumption 1 and Lem- 
ma 1, the lemma follows if 

E[A 2 (n 2 ^)]<C 
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for all sufficiently large n. Let J = B' fc A^dh , djr = diag(rA , . . . , n dk , . . . , n dus , 
. . . ,n dua )' . We write 

^ kk =j{(d^ 1 uz?\zh(d^r 1 }j > . 

We will use the inequality of Exercise 19 on page 238 of Magnus and 
Neudecker [12]. That is, 

< (jjo- 1 J{d( fc )[i m (4*\^ w )]- 1 d|*)}j / (jj / )- 1 . 

It follows that 

traced < Ai{d*)^ 1 (4 fc \4* ) )d^)}traoe{(JJ / )- 1 JJ , (JJ 1 )- 1 } 

= A 1 {d( fc )/- 1 (zf ) ,4 fc) )d(f)}trace{(JJ')- 1 }. 

Since there exists a finite constant C such that E[X 2 {d^ fc ^/~ 1 (zj , z^)dn }] < 
C for all sufficiently large n by Lemma 8, we complete the proof by showing 
that 

trace{(JJ , )~ 1 } = 0(n~ 2,ifc ). 

We write 

JJ' = n M *B' fc A fc A' fc B, + B;A< fc+ M fc+ M fc+1) A< fc+1 >B fc l {fc<a} . 

Since both matrices on the right-hand side are symmetric and positive def- 
inite, A afc JJ' > X ak [n 2d "B' k A k A' k B k ], and we have 

AittJJT 1 } < n-^AxUBlA^Bfc]- 1 } = 0(n" 2, H □ 
Lemma 10. Under Assumption 1, 
(18) p[.MX fc (H) n0Bj/o} = 0(n- 2dh i +2dk ) 

for h\ <k, k = 1, . . . , s and 



24^+2(1^ 



(19) p{MX t (H)n0B^ol=O(n- 

for fi2> k, k = 0, . . . , s — 1 . 

Proof. Since H = B<&£)B', we have X^(H) = BX^(<&d). Since *d is a 
block diagonal matrix, 

Ai(*r>) € {Aj(* fc fc)|fc = 0, . . . ,s, j = l,...,a k } 

and for Aj(<&£)) such that Aj(<&£>) = A.,- (<&&&), 

x i (* D ) = (o,...,o, x ;(* fcfc ),o,...o) / ) 



22 



W. W. CHEN AND C. M. HURVICH 



that is, the first j k _i entries are all zero. Define Em to be an event, Eht 
{A tth (* Ah ) > Ai (*«)}, < h< £ < s. We first prove (18). 

p(.MX fc (H) n®s^o} = P(X fc (* D ) + [0 Y]'), 



where the in [0 Y]' has dimension j£ x a& and Y has full rank. We have 
for h\ < k, k = 1, . . . , s, 

P(X fc ($,)/[0Y]')=p( |J E! k )< P(E c £k )=0(Y: n~ 2d ^ 

H:l<h-i ' M<hi HU<hi 

= 0(n~ 2dh i +2dk ), 

by Lemma 9. Similarly, for (19), 

p{.MX fe (H) n ®e> h2 Be ± j = P(X fc (* D ) + [Z 0]'), 

where the in [Z 0]' has dimension {q — ) x a k , and Z has full rank. We 
have for h,2 > k, k = 0, . . . , s — 1, 

P(X,($ D )^[Z0]')=P( |J £ P(^) = 0(X: n~ 2 ^) 

\:i>h 2 ' £:£>h 2 \:£>h 2 ' 



0(n 



-2d k +2d h 



8.3. Proofs for Sections 6 and 7. In this section, we will use the following 
decomposition and notation for the proofs. We write 

(20) & , AJ zz (w J -)A'& - i/Af (wj)A'b = 6'AR(w 3 )A'6 + &'AS(wj) A'6, 

where 

R(c^) = Jm^-) - *(^)/ ££ (w j )**(^) 

and 

S(o; i ) = *(^)J ee ( Wj )**(wj) " f%)- 
We will also use the following notation: 

^^'AR^jA't, 

]>>f&'AS(^)A'6, 

i=i 

m n 

^u;? d & / Af(u; i )A'&. 





An„(0 = 


1 

m n 


(21) 


MmM) = 


1 




Fm n (d) = 


1 
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8.3.1. Proof of Theorem 2. Here and in the following subsections, sup- 
pose that a £ N k , where k G {0, . . . , s}. Write d aa = d k . Note that d aa = d k . 
For 1/4 > 5 > 0, let j\f s = {d:\d- d k \ < 5}. Then for S(d) = R(d) - R(d k ), 
we have 

P(\d k -d k \>5) = P(d k e M s c n G) 



= P( inf R(d)< inf R(d)) < Pi inf S(d) < . 

Define ©i = {d : A < d < A 2 }, where A = A x when d fc < 1/2 + Ai and d fe > 
A > d k - 1/2 otherwise. Note that d - d k > -1/2 for all d G ©i. When 
> 1/2 + Ai, define ©2 = {d : Ai < d < A} and otherwise take ©2 to be 
empty. Hence, 

P(\d k -d k \>6)<p( inf S(d)<6) +P(mfS(d)<o) =o(l), 

by Lemmas 11 and 12 below. □ 

Lemma 11. Under the assumptions of Theorem 2, P(mij\fc n Q 1 S(d) < 
0)=o(l). 

Proof. Let 

U(d) = 2(d - d k ) - log{2(d - 4) + 1} 

and 

rn( j\ , G(d k ) G(d) \ 2 {d-d k ) + l ^( j \ 2 ( d ~d k ) ' 

= log — — log — — - w 2^ — 

G G(d) [ m n jT^ymr 

{-^ rn„ *\ 
5^1ogJ-(logTnn-l) \, 
m n j=1 J 

where Q = 6 / A fe ft(0)A / fc 6, as in (31), and 

G(d)=g—J2^ {d - dk) . 

n 3=1 

Then S(d) = U(d) - T(d). We have 

P( inf S(d) <0]<P( inf [7(d) < sup|T(d)| 

Following the same arguments as those on page 1635 of Robinson [17], it is 
sufficient to show that 



sup 

6i 



G{d) - G{d) 
G(d) 



o p (l). 
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Note that by Lemma 20, we have 
(22) G{d) = Cg^ d ^>C{l-e h )u, 

where £t = O p (n~ ak ). By Lemma 21, for d £ 0i 



2{d-d k ) 



\G(d)-G(d)\=£ mn (d) + M mn (d)+F mn (d)-G—Y,^ 



2(d-<4) 



tn j=l 



where £ m „, M. m „ and T mn are defined in (21). We have thus completed the 
proof. □ 

Lemma 12. Under the assumptions of Theorem 2, P(inf@ 2 S(d) < 0) = 
o(l). 



Proof. Following from the proof on pages 1638-1639 of Robinson [17], 
we write 



S(d)=log{D(d)/D(d k )}, 



where 



1 



j Mk I vv (ujj ) and v = — ^ log j . 

TThri ■ i 



Note that e v ~ m n /e. Denote 



(24) a, = < 



m r 



.2(Ai-d fc ) / ej - N 2(A 1 -d fc ) 



l<i<e", 



e < j < m n . 



By choosing A < c4 — | + ^ so that m n 1 X)j=i( a j — 1) > 1 f° r a h sufficiently 
large m n , we have 



J 7 



> 1 
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Now, by (20), 



1 



— EK-i) 



m 



n i=i 
(25) 



I vv (ujj) 



n j=i 



I vv (Uj) I vv (iOj) 

guj -2d k ~ VM {uj^)A'b 



^b'AR(ojj)A'b 
+ m n j^p 1} 6'Af(^)A'6 

b'ASju^A'b 
+ m np[ j ] VAf{ujf)A'b' 

We will show that all three terms in (25) are o p (l). For the first term, we 
begin by showing that 



(26) I vv {uj) = b'AR{uj)A'b + 6'AS(w i ) A'b + b'Ai{uj)A!b = O p (, 



-2d,, 



Let ¥Lhe(ujj) denote the (h,£)th block of H(ujj). By Lemmas 16 and 18, 



(27) 





' O p (n 2d "- dh - di ujJ dh ' 




h,£<k, 


b'A h K h i(ujj)A' e b = < 


o p ^~ d ^r pl2 \ 




h,£>k, 




O p {n d ^ujj dh - dl j- 




h<k,£>k 


= O p {u>- 2d *j- p/2 ). 







Also, by Lemmas 16 and 19, 
b'AS(uj)A'b 



k-l 



~ d t n d-2k-dh- 



\h,e=o 

O p {u] 2dk {j 2dk - 2dk -^ + 1 + 



+ E-7" + E -7 

h<k,£>k 



~ d i n dk-dh 



By (54) in the proof of Lemma 20, 6'Af (uj)A'b = O p {ujJ 2dk ). Thus, the 
bound in (26) follows. Together with Lemma 20, we have 



I vv (u)j) 



o P (i) 



and 
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Thus, the first term of (25) is 



i m " / Quj- 

— y>i-i) (1-17 3 

■Tin pi J \ 



(28) 



2d k 



I vv (iUj ) 

b'Ai{u;.)A>bJ~^^ 



1/2 



since X)j=i Q j = 0{m^ dk A ^ + mlogm), by Equation 3.24 of Robinson [17]. 
Applying (27) and (26), the second term of (25) is 



1/2 



O, 



-p/2 



o P (l), 



n i=l 



V/=l 



by the same argument as for (28). The third term of (25) is bounded by 



m n ^ aj > VAi{u-.)A>b 

,7-1 ** 



+ 



??? 



1 ^ . iN 6'AS(u; 7 )A'& 

K-i)t7 

" j=[e«]+l 



VAffa)A'b 



Following from (24) and the lower bound of b'Af (u)j)A'b in Lemma 20, the 
first term of the above equation is 



« j=l 



by (ii) of Lemma 21, because > A — dk> — 1/2. We will complete the proof 
by showing that 



(29) 



1 



E 



2(A 1 -d fc ) & 'AS(a; i )A'6 
6 / Af(wj)A'6 



o„(l). 



Note that e 1 ' ~ m n /e. Following the similar computation for (55), 



E 



1 



E 

" J=[e«]+1 



2(A 1 -d fc )g 



q| n -4d k +2d h +2d em 2(2d k ~2A 1 ~l) j4Ai-2d h -2<i £ 

\ j=Sm n 
0{uj 4d 2d h -2d lm -ly 
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Hence, 



1 



E 



j=[e«]+l 

By Lemma 16, we have 
1 



P « fe " dh ^^n 1/2 )- 



V m » i= [e-]+l Vm ^ y 



O v (m 2 n k ~ dh ~ dl ~ Xt2 ) 



h,£<k, 

O p (u 2 Jy d ^mn 1/2 ), h,£>k, 

p p (cv^mn 1/2+dk - dh ), h<k,£>k, 
= o p (l). 

Equation (29) follows from the triangle inequality. □ 



(30) 



8.3.2. Proof of Theorem 3. By Theorem 2, 4 satisfies 
8R(d k ) dR(d k ) d 2 R(d) 







dd dd 
where \d — dk\ < \dk — 4| ■ Let 

m„ 

-1/2 



+ 



dd 2 



-(4 -4), 



2m n 1 ' 2 J2 Vj^eeiUj) ~ £), I/j = log j ^ log j 



and 

where 
(31) 



3„ = ^6'A fc *t , (0)Z n *t(0)A' fc 6, 



6? = &'A fe f+(0)A;6 = 6'A fe *i'(0)S*l (0)A' fc 6 



and isagx a fc submatrix of *t(u;) = [*q(w) • • • in ( 8 )- We 

show in Lemmas 13 and 14 that 



(32) 
and 
(33) 



d 2 R(d) P 



dd 2 

l/2 dR{d k ) 
dd 



3n + O p (l). 
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From Lemmas and 8 of Hurvich and Chen [8], the (u, v)th entry of Z n 
satisfies 

Using a similar computation for the variance above and Eq. (53) in the proof 
of Lemma 19, we obtain 

Using the Cramer- Wold device, we have vec Z n —> vec Z ~iV(0, 4$ p S (g) S). 
By Lemma 15, b^>b, thus, 

(34) 3,4 ^f(0)Z^l(0)A;& ;= 

b'A k ^ (0)^1(0) A' k b 

Let tp = (<p 1 ,...,<p q )' = ¥%(0)A' k b. Then 3|6 is a normal random variable 
with mean zero and variance 

Thus, 3 is independent of b and 3 ~iV(0, 4<J> P ). Together with (30) and (32)- 
(34), we have proved the theorem. □ 

Lemma 13. Let d be such that \d — d&| < |<ifc — T/ien under the 
assumptions of Theorem 2, 

d 2 R{d) v , 

Proof. Define 

i m n 

G a (d) = — £(logw 5 )°wf J w ( Wi ) 

and 



F a (d) = — £>g j) a a;f / TO (^), E a (d) = — £>g j) a j M / w (^). 
n j=i n j=i 

Then 

d 2 R(d) _ 4{G 2 (d)G(d) - G\(d)} _ 4{F 2 (d)F (d) - F?(d)} 
dd? &{d) F${d) 

4{E 2 (d)E (d) - E\{d)} 
El{d) 
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We first show that 

(36) F a (d) = F a (d k )+o p (l), a =1,2,3, 

by showing that E a (d) = E a (d k ) + o p (n 2dk ) for a = 0,1,2. Let M = {d : 
log 3 m n x \d — d k \ < e}, where e > is fixed to be such that 2e < log 2 m n 
with a proper n. Following the same line of proof as on page 1642 of Robin- 
son [17], for 7] > 0, 



P[\E a (d)-E a (d k )\> V ( — ^ 



(37) 



V n 



<P[G(d k ) > 2^(logm„) 2 - a J + P(log i m n \d-d k \>e). 
The first probability is bounded by 

p(\G(d k ) -g\> ^L(logm n f- a ) +p(g> ^-{\ogm n f- a 



4ee 



4ee 



Both probabilities in the above equation tend to for e sufficiently small 
since \G(d k ) — G\ = o p (l) and Q < C, by Lemma 20. To show that the second 
probability in (37) tends to 0, we only have to verify that 



sup 

eirW,, 



G(d) - G(d) 



G(d) 



o p (log 6 m n ). 



From (22) and (23) in the proof of Lemma 11, 



sup 

0irW s 



G(d) - G{d) 



G(d) 



< sup 
6i 



G(d) - G(d) 



o p (m n e ). 



We have established (36). Combining this with (35), we have 
d 2 R{d) 4{F 2 (d k )F (d k ) - F 2 {d k )} 



3d 2 

By Lemma 21, 



F a (d k 



g—J2iog a j 



m 



£iog a j 



G(d) 
th 

+ o p (l) as n —> oo. 

I vv (Uj) 



U)~. 

J 



-2d h 



< \og a m n \C mn (d k ) + M mn {d k ) + T mn {d k ) -Q\ 
= O p {m~ £ log" m n ). 
By the same reasoning as that used in (4.10) of Robinson [17], we obtain 



8 2 R{d) 
3d 2 



m 



^log 2 j 



m, 



glogj (l + 0p (l)) +0p (l)^4. □ 
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Lemma 14. Under the assumptions of Theorem 3, 

m y2^)_i 6 / Afe *t* ( o )Zre *t (0)A;6 = 0p(1)) 



where 



Z n = 2m; 1/2 ^j/ J (/ ee (w i ) - S) and Vj = log j ^logj. 

-i TThn ■ -i 



Proof. Note that 



&R(d) _ 2 g Vjl m (u}j) 



3d m n j± x ioj 2d G(d) ' 
Since G(d k ) — G = o p (l), by (23), and X)j=i u j = 0> we have 

ll/2 aflW0 = 2m -i /2 g ^ / ^-M _ A (1 + 0p(1)) 



3d n ^ v Vgo;-7 z ' 



and 

m l/ 2 ^)_l 6 / Afe ^ (0)Zn ^t (0)A / 6 



(38) 

+ 2m" 1 / 2 1 J] ^(u;f fc 6'A fc S fcfc (^)A' fc 6) - &'A fe *t*(0)Z n *t (0)A j 

+ Op(l), 

where S«(tUj) is the (/i,£)th block of S(cjj) defined in (20). Let 

Mt! ] (d) = — f>f fr'A^S^A^. 

m n . =l 

The first term of (38) is then 

i=i Gu< k { J 



u,v — 
v^k 



2ml/ 2 log m„ 



G 



Cm n (d k )+F mn (d k )-G+ £ ^tfe?(d* 



u,i> = 
vj^k 
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by (20), (21), Lemma 22 and Assumption 3B. Since 

s tt ( Wi ) = ii - e ^r Mfc *i*K)[i(^) - 

= W? - a **f (0)[/( Wi ) - S]*I(0) + O p (u,-7 2d * + 0, 
the second term in (38) is 

by Lemma 16 and Assumption 3B. We have shown that both terms on the 
right-hand side of (38) are o p (l) and, hence, have completed the proof. □ 

Lemma 15. Under Assumption 1, the matrix X = X.{I m (y t ,yt)} satis- 
fies 

X 5. X = B'X{K(UU' + VV%K'}, 

where X{K(UU' + VV') B K'} is the matrix of normalized eigenvectors of 
K(UU' + VV')oK' in Lemma 6 and U,V are defined as in Lemma 1. 
Thus, X is a continuous function with respect to vec(U, V). 

Proof. It suffices to show that 

(39) \\X k (H)-X k \\ F = O p (n-^) 
and that the eigenvectors of H satisfy 

(40) Xj(H )^4(vec(U,V)), j = l,...,q, 

where are continuous functions of vec(U, V). 

We first show (39). Note that since both H and L m are symmetric, we can 
assume that Xj Aj(H) > 0. We have 

||X fe (H) -Xfc||| < a k max||xj(H) - Xj\\ 2 < Cmaxsin 2 6»(xj,Xj(H)) 

<C||sine{^(X fc (H)),^(X fc )}|| 2 ,, 

by the definition of the sin O bound. Equation (39) follows from Lemma 4. 

Next, we derive (40). Since Xj(H) = ~B'xj(*&D), it is sufficient to show 
that 

Xj(*2j) C,(vec(U, V)), j = l,...,q, 
where <6 are continuous functions of vec(U, V). Let = d„ €?£)d~ . Then 

* D = X'(* D )A(* fl )X(* D ). 
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First, note that the eigenvalues of are distinct, with probability 1, by 
Okamoto [14]. Since both $>d and are block diagonal matrices, we have 

* fl = d- 1 X'(* D )A(* D )X(* D )d- 1 = X'(* D )d^ 1 A(* I) )d- 1 X(* D ). 

This implies that 

X'(* D ) = X'(*s) and A(* D ) = d^ 1 A($ D )d^ 1 . 

We now let K be defined as in Lemma 6 and rewrite U n and V n in (10) as 

U n = [U' o ... U' ]', V n = [V' o ... V; 



n,sl ! 



where \J n ,k and ~V n ,k are a k x m matrices. Since K is a block diagonal 
matrix, we have 

$ D = K' diag(U n , U; i0 + U n , U; i0 , . . . , U n , s U^ + U^U^ J K - 

It follows that 

Xj (* D ) = : = ff(vec(U n , V n )) 9(vec(U, V)), 

where d-(-) is a continuous function because the eigenvalues of 4>£) are dis- 
tinct with probability 1 and 

Xj (H) = B / Xj (*D) -^B^(vec(U,V)) = 4(vec(U,V)). □ 

Remark 3. Let X = [X ... XJ. Since X{K(UU' + VV%K'} is a 
block diagonal matrix, ± k = B'Xj,{K(UU' + W%K'} G B fc . 

8.3.3. Proof of Theorem 4. In case of no cointegration, we have 
d = Ai(AA') > ||fe'A|| 2 = (b'AA'b) > \ q (AA') > C 2 

and 

Ci = Ai(AA')Ai(f t (0)) > g = b'Af\o)A'b > \ q (AA')\ q (f\o)) = C u 

where C±, C 2 , C\ and C 2 are positive constants. Furthermore, by Assump- 
tion 2, 

b'Ai (uj)A'b - uj 2d0 g = O p (u;J 2do+p ), 

for 1 < j < m n . Following along the lines of the proofs of Theorems 2 and 3, 

we have ml/ 2 (d a a — d aa ) N(0, 3>p/4) for a = 1, . . . ,q. The theorem follows. 
□ 
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9. Technical lemmas. We will need the following two lemmas: 

Lemma 16. Ifb = Xa, wherea£N k , then, under Assumption 1, 

b'A h = O p (n- d » +dk ) 

for h < k, k = 1, . . . , s and 

b'A k = O p (l), 

for k = 0,...,s. 

Lemma 17. Ifb = Xa, where a£N k , then, under Assumption 1, 

\\b'A k \\>C(l-e k ), 
where C > and e k = O p (n~ ak ), k = 0,...,s. 

PROOF of Lemma 16. Since X(H) is an orthogonal matrix and .MX(H) = 
R 9 , we have 

s 

(41) 6 = £x*(H)q, 

1=0 

where 

c i = X' i (R)b = O p (n-^- d ^), 
by Lemma 7. Furthermore, for I > h, 

E[||XKH)A /t ||]=E[||X^(H)A,||l { ^ x ,(H)ce J> ^ J }] 

+ E[||X^(H)A /l ||l{_ MX( ,(H)n® j < h B j ^o}] 

< + E[||X^(H)A h ||l{_ MXf (H)n© j < h £? J ^o}] 

= E[trace 1 / 2 (A / h X,(H)XKH)A h )l { ^ x ,(H)nffi J <^^o}] 

(42) 

< E[trace 1 /2( A / /l A fe ) trace 1 / 2 (X,(H)X^H))l {A1XK H)ne J<h B^o}] 



1/2 , , . , 

a/ \\A h \ 



1/2 

pjjKx,(H)n0e^o|' 

j<h 



= 0{n~ dh+di ), 
by Lemma 10. For £<h, 

(43) E[||Xj(H)A fc ||]=0(l). 
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We have, for h < k, 
b'A h = J24X.' e (H)A h 

£=Q 

= ]T c£Xj(H)A fc + ]T 4Xj(H)A fc 



l:£<h l:£>h 

-d e +d k + 



Op E 



n -d h +d t -d t +d k _|_ n - 
£:h<£<k £:£>k 



d h +d t -d k +d t 



= O p {n- dh+dk ). 

For h = k, the above equation is of O p (l) since = andE[||X' fc (H)Afc| 
0(1). □ 

Proof of Lemma 17. Note that 



(44) 



4X' fc (H)A fe + ^ 4X',(H)A fe 

£=0,i^k 



> 



||4X' fc (H)A fe | 



Using (41), we have 



i = n&n 2 =]Tiix,(H) Q |i 2 =]r 



£=0 

and 
(45) 

by Lemma 4. Thus, 
(46) 

By (45), (42) and (43), 



E 



2 = O p (n- 2 ^), 



l-O p (n- 2a *). 



£ c£x£(H)A, 



+ E 

£=0,£^k 



(47) 



1/2 



< £ Nil 2 E l|Xj(H)Afc| 
= O p (n- a fc). 



Furthermore, if A^Xfc(H) = B)~, then there exists an x orthogonal 
matrix D such that 



X fe (H)=B fc D 
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since both X&(H) and are matrices with orthonormal columns. We have 
||c fc || 2 = trace{ C , fc DB , fc A A: (B / fc A fc )- 1 (A , fc B fc )- 1 A' fc B fc D / c fc } 
< ]|(B' fc A fe )- 1 || 2 ||4DB^A A; [| 2 
= ||(B' fc A,)" 1 || 2 ||4X' fc (H)A fc || 2 . 

It follows that 

H4x' fc (H)A fc || 2 > ikb'.a,)- 1 !!- 2 !!^! 2 = c(i - s k ), 

where 5 k = O p (n~ 2ak ), by (46). By (44), (47) and the above equation, ||6'A fc || > 
C{\ — 5k — ik) , where ik = O p (n~ ak ) . We have thus completed the proof. □ 

Lemma 18. Let R a b(u>j) be the (a,b)th entry ofH(ujj), 
®\Rab(uj)\ < C\l- e -^l\-( d ^+ d ^j-p/\ 

a, b = 1, . . . , q and 1 < j < [n/2] 

under Assumption 2. 

Proof. Let J Za {oJj) be the jth element of J z (ujj), the discrete Fourier 
transform of Zt- By (4), 

(48) JzM)=ib j z ab (^ 



b=l 



where 



/2ttE|^ 

Hence, 



i n / oo \ 



R ab (Uj ) = J Za {ujj ) J Zb (Uj ) - J2 *«* ( W i ) J eu {"j ) ^bv (WJ ) J £v (Wj ) 

(49) 

(*au(w--)*6w(^)(-4.au,i^6« ! j ~ B u jB v j)), 



q 

E 

u,v=l 



where 

(50) ^ = t°"M and B Utj = J £u (uj). 

From Lemmas 9 and 10 of Hurvich et al. [9] , 

E\A auJ -B uJ \ 2i < c( r p^-l^D^Uj-utfdw)' 

\J—K ^au\^) / 



(51) 



36 W. W. CHEN AND C. M. HURVICH 

By the Cauchy-Schwarz inequality, 

"I - Buj^Afoj B Vt j) + B v j(A au j B u j)\ 
< 3(E\A auJ - B uJ \ A E\A bvJ - B vd \ A ) 1 / 2 

(52) 

+ {E\B u ^\ 4 E\A bV)j - B^) 1 ' 2 
+ (E\B v j\ A E\A auJ - B ud \y 2 

<c[(r 2p r 2p ) 1/2 + (r 2p ) 1/2 ] = cj^. 

We have, from (49) and Assumption 2, 

<? 

E|i? a6 |< ]T VauiUjWbviu^AaujA^j - B UJ B V J 2 Y' 2 
u,v=l 

u,v=l 

<C\1 — e~ iu) 'i \~( daa+d t>b)j-p/ 2 
where the constant C does not depend on n. □ 

Lemma 19. Let S ab {ijo) be the (a,b)th entry ofS(u>j). Then for 1 < j, k < 
[n/2], 

r C |(l- e -^)(l-e-^)|-(^+^), \j-k\<p, 

E\Sab(wj)Sab(Wk)\ < \ . 

(C/n, otherwise, 

under Assumptions 1 and 2. 

Proof. Note that EJ ee (wj) = £ and 5 a6 (w i ) = Y? u ,v=i *au(^)* 6 „(6^) 

X (^ee ,m>( w ,j) ~~ °~uv)- Now, 

Q 

E|5 a6 (a;j)5 a b(a;fe)| = ^ *aui(^)* a u 2 (^*:)*6t;i(wj)*6« 2 (^fc) 

«l,«2,"l,f2=l 

X E[(/ eejUlt , 1 (wj) — 0" U1U1 )(J ££)U21 , 2 (wfe) — cr U2t , 2 )]. 

Note that E(J £ii ( Wj ) J e „(w fc )) =0, 1 <j,k < n/2, and E( J £u {u: j )J £v (w fc )) = 
if \j — k\ > p and 

(53) E(J £u ( Wi )^K)) = ^(-ir fc ( p l{U- fe |<P}> 
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where 

see [9]. Hence, 

^[(-^ee,«iUl \y)j) ®U\Vi)(JsE,U2'V2v^k) '-'"'"2^2)] 
= E[I ££jUlVl {uJj)I e£ ^ U2V2 (u) k )\ 0~ UlVl O~ U2 v2 

= cum( J £ui (uj), J £u2 (u k ),J £vi (uj),J Ev2 (a;*,)) 
+ E(J £ui {u 3 )J £v2 ( Wfc ))E( J £U2 {u>i)J evi (u k )) 
= Cl{\j-k\< P } 

because (53) and the cumulant is under Assumption 1. We have, by (5), 

E\Sab(Vj)Sab(Vk)\ 

g 

<C ]T (1^(^)11^2(^)11^(^)11^2(^)1) L {\j-k\<p} 

Ml,«2,«l,V2=l 

<C|1 -e-^l-^+^^li _ e -K|-(^+^) 1{|j _ fc|<p} . n 

Lemma 20. Under Assumptions 1 and 2, 

b'Af (uj)A'b - b'A k f kk (u>j)A' k b = O p (LjJ 2d "(j d "- d ^ + uf' dk+1 )) 

and 

b'A k f kk (u; 3 )A k b - Gu] 2dk = O p (u;- 2dk+p ) 
for l<j< m n . Furthermore, there exists a constant C such that 

b'Af (uj)A'b > Cuj] 2dk (l - e k ) 
for l<j< m n and two constants, C\ and C2 , such that 

Ci>G>C 2 (l-e k ), 

where e k = O p (n~ ak ). 

Proof. Since, by Lemma 16, 

f O p (ujj 2dk j 2dk - d ^), h<k,£<k, 

< O p (uJ dh - dl ), h,£>k, 

O p (w: 4 " d 'i^- rf "), h<k,£>k, 



(54) b'A h i M (io 3 )A' e b-- 
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we have 

s s 

h=o e=o 

= b'A^ki^Kb + OpiuJ 2 ^^-^ + u~ 2dk+1 + u- dk - dk+1 ) 

= 6'A fe f fefe (^)A^ + O p (a;- 2d H/ fe - dfc - 1 +^' ! ~ dfc+1 ))- 

Since (7) and Assumption 2 imply that f kk (uj) =f^(0)w- Mfc +O{oj~ 2dk+P ) 
as uj — > 0, we have, by Lemma 16, 

b'Af kk (uj)A k b' = b'A k fl k (0)A' k bu;J 2dk + O p (\\b' A k \\ 2 ooJ 2dk+p ) 

= guJ 2dh +O p (coJ 2dk+p ). 

We have shown the first two equations of the lemma. For the third equation, 
we have, by (54), 

b'Af(uj j )A'b = b'( ]T A h f M (u j )A'^\b + O p (u;J dk - dk+1 ). 

\h4=0 J 
By Assumption 2 and Lemmas 16 and 17, 

b'(j2 A h f M (u j )A')ib = b , (Y^ ^^A fe 4(0)A^6 

\h,e=o / \h,e=o / 



+ O p [uj dh+d * +p £ b'A h A! t b 



k 



hl=0 



k 

>uj 2dk \ min {f\0)} E b'A h A' e bA h + O p (ujj) 
hJ=0 



>Ca;7 Mfc (l-e fc )+O p (^). 

For the last inequality, 

G = b>A k fi k (0)A' k bujJ 2dk > W 7 Mfc A min (4(0))||6'A,|| 2 > Cuj 2dk {l - e k ), 
by Lemma 17. The upper bound for Q is due to the fact that 

a<A max (4(0))||A fc ||||6||=A max (4(0))||A fc ||. □ 

Lemma 21. Let C mn (d), M mn (d) and J- mn {d) be defined as in (21). 
Then if d — d k > — \, there exists an e > such that 

(i) C mn (d) = o P (u 2 J- 2dk m-% 
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(ii) 



M m M=oMt 2dk ^n £ ), 



(iii) T mn {d) - Q^~ 2d * = o p (u;^m-% 

under Assumptions 1 and 2. 

Proof. We will only prove (ii); (i) and (iii) can be shown in a similar 
fashion using Lemmas 16, 18 and 20. Let S^(d) be the (h,£)th block matrix 
of S(d). By Lemma 19, 

2 

E 



1 m„ p+j 

^ 2d—d h —df 2d—d h —de 



"In P+j 



1 f/tn ^ i j 



4d-2d h -2d£ 



'0(n Mh+M «- 4d m- 2 logm n ), 
0«T 



Ad-2d h - 2d e < -1, 
Ad-2d h - 2d t > -1. 



Hence, we have 



(55) 



Let 



' O p (n dh+d «- 2d m- 1 log 1 / 2 
O p «-*-^m^ 1/2 ), 



2d-d h -d e < -1/2, 
2d- d h -dt> -1/2. 



I 1 . . \ 



By Lemma 19 and (55), we have, for h,£ < k, 
M&!\d) 

(56) { O p (^-^m- l - 2d + 2d " log 1 / 2 m ra ) 



.°p( w 4t 



2d-2d,.^ n - 1 /2+e 



■m 



), 



2d-d h -d £ < -1/2, 
2d-d h -d e > -1/2, 
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where e > 0. By the same lemma and (55), we have, for h,£>k, 

(57) Mt'f 1 (d) = O^^m- 1 / 2 ) = O p (^;Xr 4 "^n V2 ) 
and for h < k, £> k, 

(58) (opiu^-^m-^+^log^mn), 2d- d h -d e < -1/2, 



o p (ult 2dk ^- de ^n 1/2+dk - dh ), 2d-d h -d e > -1/2. 



Hence, 



h=oe=o 

since 2d k - d h - d e > in (57) and -1 - 2d + 2d k < in (56) and (58). □ 
Lemma 22. Under the assumptions of Theorem 2, if d — d k > —-g, then 

and 

'O p (uj%- 2d *m- 1/2 ), h = £ = k, 



_o p (uJ 2d n 2dk m Tl l l 2 e ), otherwise, 



Furthermore, under the assumptions of Theorem 2, 

T mn {d k ) ~Q = O p (m d n k - dk -' + ufc dh+1 ), 

where the O p (mn k dfe_1 ) term is vacuous if k = and the O p (t0m n a!fc+1 ) term 
is vacuous if k = s. 

Proof. This lemma is a corollary of Lemma 21. □ 
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